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We present first results from a number of experiments conducted on a 0.53 kg cylindrical dumbbell-shaped 
sapphire crystal. This is the first reported optomechanical experiment of this nature utilising a novel modifica¬ 
tion to the typical cylindrical architecture. Mechanical motion of the crystal structure alters the dimensions of 
the crystal, and the induced strain changes the permittivity. These two effects result in parametric frequency 
modulation of resonant microwave whispering gallery modes that are simultaneously excited within the crys¬ 
tal. A novel low-noise microwave readout system is implemented allowing extremely low noise measurements 
of this frequency modulation near our modes of interest, having a phase noise floor of -165 dBc/Hz at 100 
kHz. Fine-tuning of the crystaUs suspension has allowed for the optimisation of mechanical quality factors in 
preparation for cryogenic experiments, with a value of Q = 8 x 10 7 achieved at 127 kHz. This results in a 
Q x / product of 10 13 ; equivalent to the best measured values in a macroscopic sapphire mechanical system. 
Results are presented that demonstrate the excitation of mechanical modes via radiation pressure force, allow¬ 
ing an experimental method of determining the transducer’s displacement sensitivity df/dx, and calibrating 
the system. Finally, we demonstrate parametric back-action phenomenon within the system. These are all 
important steps towards the overall goal of the experiment; to cool a macroscopic device to the quantum ground 
state at millikelvin temperatures. 


I. INTRODUCTION 


T he field of optomechanics is producing many excit¬ 
ing results, from extremely precise sensors with ap¬ 
plications in a wide variety of fields, including sensitive 
detection of previously immeasurable signals 1 ^ 7 , quan¬ 
tum information processing 8,9 ,and tests of fundamental 
quantum theories, including potential tests of quantum 
gravity 10,11 . The majority of work in this area is based 
around microscopic or mesoscopic resonators, which have 
the advantage of extremely low mass and therefore rel¬ 
atively large optomechanical coupling factors, allowing 
for a plethora of interesting physics to be investigated. 
In particular, the quantum regime of mechanical motion 
(the so-called “standard quantum limit”) is far more ac¬ 
cessible for these extremely small resonators. Great chal¬ 
lenges arise when dealing with macroscopic resonators, 
and it has been shown that state-of-the-art technologies 
are required to approach observations of the quantum 
world in such resonators, due to their significantly larger 

19 

masses ’ . 

The essential requirement for the observation of op¬ 
tomechanical effects is to generate coupling between 
phonons and photons within a resonant structure. This 
coupling can take many forms, from variations in path 
length of a Fabry-Perot resonator via a moveable mir¬ 
ror (such as the LIGO experiment) 11,14 , oscillations of 
a cantilever capacitively coupled to two electrodes 15 , 
“membrane-in-the-middle” type experiments 16 , or as in 
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the present work, crystalline resonators, which can sup¬ 
port both optical and microwave resonant fields and me¬ 
chanical modes 17,18 . All of these systems function via 
mechanical effects causing a disturbance to an optical 
field, which can be measured as a frequency shift. In the 
case of a macroscopic whispering gallery mode (WGM) 
resonator, its frequency is sensitive to changes in path 
length (i.e. the circumference and height of the cylin¬ 
der), and strain in the crystal lattice, both of which are 
induced by mechanical motion. The latter effect is in 
fact the more dominant mechanism 12,13 . These lattice 
vibrations can be caused by thermal phonons within the 
resonator at room temperature. 


II. SAPPHIRE RESONATOR: GEOMETRY AND 
MODES OF VIBRATION 

The sapphire resonator studied takes the form of a 
dumbbell-shaped cylinder and will herein be referred to 
as the “Split Bar” (SB) resonator, with its dimensions 
outlined in Fig. la. It is suspended via a wire loop 
around the central “neck”. This shape was chosen to 
isolate the point of suspension from the mechanical 
motion in order to maximise mechanical Q-factors. The 
crystal is high quality HEMEX grade sapphire and 
was grown using the heat exchange method by Crystal 
Systems, USA, cut to dimensions and then optically 
polished. 

The SB is both an electromagnetic and mechanical 
resonator. Its cylindrical shape permits high Q-factor 
WGMs to be excited at microwave frequencies, and via 
the mechanism explained above the coupling between 
acoustic and electromagnetic modes required for the 
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FIG. 1: (a) Graphical representation of the SB 
resonator modelled using FEM. Lengths shown are in 
millimetres. Normalised displacement field of the 95 
kHz (b) and 127 kHz (c) resonant mechanical modes. 


observation of optomechanical effects is achieved. 

The SB resonator was modelled using Finite Element 
Modelling (FEM) software (COMSOL™) to identify 
suitable mechanical modes within the crystal. Given the 
non-uniformity of the dimensions, the eigenfrequencies 
of the resonant modes were solved for varying neck di¬ 
ameters, ranging from the limiting cases of two separate 
cylinders of height 27.4 mm, to a single uniform cylin¬ 
der of height 63.75 mm, in order to identify the correct 
families of the modes. 
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FIG. 2: Three relevant mechanical mode families 
resonant frequencies as a function of neck diameter. 
The SB resonator has a neck diameter of 14.92 mm. By 
examining the two limiting cases of the SB (i.e neck 
diameter = 0 mm; corresponding to two separate 
cylinders, and neck diameter = 27.4 mm; corresponding 
to one large cylinder, the origin of each mode can be 
correctly identified.) 



From Fig. 2, we can see that the in-phase and out of 
phase breathing modes (127 kHz and 125 kHz modes, 
respectively) are degenerate for the case of two separate 
crystals, but when they become joined, their frequencies 
tune in opposite directions. These modes become very 
different in the limiting case of a full cylinder. Each one 
of the modes in Fig. 2 has been detected with frequen¬ 
cies in very good agreement with the modelled values 
(< 1%). However, this article will mainly deal with the 
95 kHz mode (black curve) and the 127 kHz mode (red 
curve), as they have the two highest Q-factors. The low 
Q-factor observed for the 125 kHz mode is predicted by 
the modelling due to a large amount of displacement at 
the suspension point on the neck of the dumbbell, result¬ 
ing in larger suspension losses. 

From the FEM, it is possible to obtain deformation 
gradients (^, , ^) and total displacement curves for 

each of the modes, allowing an estimate of the WGM 
frequency sensitivity to each different mechanical mode 
family. For simplicity, this is done using a model of only 
one cylindrical end of the SB, at the mode of interest’s 
corresponding eigenfrequency in this limiting case, as 
depicted in Fig. 2. 

To determine the magnitude of frequency variations due 
to mechanical vibrations (^) of any particular WGM, 
we expand on the method used by Locke et al. 12 . The 
frequency of WGMs in cylindrical sapphire, assuming the 
crystal axis (c-axis) is aligned with the z-axis, is depen¬ 
dent upon four variables; the permittivities of the crys¬ 
tal perpendicular and parallel to its c-axis, e± and £y, 
respectively (sapphire is an anisotropic crystal), and the 
crystal’s dimensions (diameter, D and length, L ). 

The frequency sensitivity of the WGM can be calcu¬ 
lated from 


jf z = V (( M rPer + M+PeJ K e± +p D ) 

- M z Pe z K e]] - PL, (1) 

where Pi are the normalised tuning coefficients of the 
WGM frequency with respect to the variable i (i.e. pi = 

j), v is Poisson’s ratio, K e . represents the strain 

dependence of permittivity (i.e. K e ^ = and K e± = 

2 F 7 )’ an d Mi is the displacement modification factor, 
and represents the overlap of electric field (Ej(r, z, <ft)) 
and strain (S^(r, z, </>)) in the i direction; 
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( 2 ) 


where Ni is a normalisation constant. 


Frequency sensitivity to changes in permittivity only 
occurs in the electrical component of the resonant 
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microwave mode. This is due to the dependence of the 
E-field on 6 compared to a B-field’s dependence on /i, 
a direct result of Maxwell’s equations. WGMs have 
two main polarisations - WGH modes with dominant 
E z , H r and field components, and WGE modes 
with dominant H Zl E r and E^ components. Therefore, 
transduction in the split bar will be primarily due to the 
strain induced in the direction of the WGM’s E-fields. 
Since v = 0.3 for sapphire, WGH modes are more 
strongly affected by mechanical motion. 
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FIG. 3: Overlap of WGH 1511 E z field and deformation 
gradient used to calculate M z modification factors in 
equation ( 1 ). 


Figure 3 depicts the method by which the M z factors 
are calculated for the 95 kHz and 127 kHz mechanical 
modes. The overlap of the electromagnetic and strain 
fields determines the magnitude of WGM frequency vari¬ 
ations. It should be noted that the 95 kHz mode’s de¬ 
formation gradient is an antisymmetric function about 
Z = 0, and hence the mean strain is zero, compared to 
the symmetric 127 kHz mode. This results in effectively 
zero strain contribution to ^ in equation (1) for the 95 
kHz mode, hence the main component of transductance 
for this mode is expected to be a result from the physical 
change of boundary conditions. 

Following the treatment used by Tobar and Mann 19,20 , 
a mode matching technique is used to determine the pi 
values and the E-field distributions in each of the cylin¬ 
drical coordinate directions. This can be done for both 
WGH and WGE modes with any number of azimuthal 
maxima. Then, taking the deformation gradients for the 
two mechanical modes produced in FEM, Mi values are 
determined to finally produce an estimate for as shown 
in Table I. 


Ho/ 27r (kHz) 

M z 

M r 


\df/dz\ (MHz/pm) 

94.97 

0* 

0* 

0* 

0.085 

127.07 

0.38 

0.0017 

0.024 

0.18 


TABLE I: Calculated values for the WGHi 5 ? i 5 i mode 
and the 95 kHz and 127 kHz mechanical modes. * a 
value of zero is due to the asymmetry of the 95 kHz 
mode strain curves. 


III. EXCITING MECHANICAL MODES AND 
MEASURING Q-FACTORS 


Initially, a piezoelectric actuator was used to excite the 
SB and therefore locate the resonant mechanical modes. 
The piezo was mounted to the SB’s copper housing 
on a lid at one of the end faces, and is not in direct 
contact with the crystal. This setup is shown in Fig. 4. 
The piezo driving frequency was chirped over a small 
frequency range (on the order of 1 kHz) to conduct 
a relatively broadband search for resonances around 
the predicted frequencies. However, under vacuum the 
only path for the acoustic excitation of this form is 
through the suspension wire and if the length of this 
wire is optimally tuned such that the driving frequency 
is non-resonant with the wire, the transmission of energy 
to the bar will be minimal and the modes difficult to 
excite. In addition to this, driving the piezo with a 
periodic chirp results in the reduced contrast of the 
resonant peak viewed on a network analyser due to 
non-resonant excitation caused by shaking of the copper 
housing. 

An alternative method for mechanical excitation is 
applying radiation pressure force (RPF) to the crystal 
by modulating the power of the AM microwave signal at 
the mechanical resonance frequency (see Fig. 10). Once 
the AM is switched off the SB will continue to resonate 
mechanically with its amplitude ringing down with a 
characteristic time constant, r, which can be used to 
calculate the Q-factor of the mode (Q m = tt rf mec h). 
The results of this ringdown technique are shown in Fig. 
5. 

When used together, these two techniques provide an 
effective and simple method for the excitation of me¬ 
chanical modes and Q-factor measurements - the piezo 
is chirped to find the exact location of the resonance 
(within one bandwidth, on the order of mHz), and the 
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FIG. 4: A diagram of the SB suspended inside a copper 
cavity. The length of wire from the point of contact with the 
crystal neck to the position of the niobium clamp, labelled 
in blue, is the length that must be tuned. 
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FIG. 5: Mechanical Q-factor can be obtained from 
fitting the ringdown of the mechanical mode, measured 
by the output of the frequency discriminator. 
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FIG. 6: Measurements of maximum peak amplitude of 
readout system output (proportional to displacement) 
to varying amounts of RPF. 


input microwave signal is then AM-modulated at this 
resonant frequency for a period of time longer than the 
ring down time of the mode (to allow full excitation). 
The AM-modulation is then switched off, and finally the 
amplitude of the resulting peak is tracked over a time 
period of about 200 s. Fig. 5 shows a measured value of 
r = 1/0.0052 s for the 127 kHz mode, corresponding to 
a mechanical Q-factor of 7.7 x 10 7 . 

By varying the amount of modulation on the incident 
power and measuring the peak amplitude immediately 
after AM-modulation has been switched off (provided 
RPF has been applied for a sufficiently long period of 
time that the resonator is in a quasi-static regime), one 
can experimentally determine a value for df /dx for each 
of the mechanical modes, and compare with predicted 
values, as will be shown. 

Treating the mechanical sapphire resonator as a stan¬ 
dard harmonic oscillator we have: 

d 2 x(t) _ dx(t) o / x _ , , 

ra e ff 2 +rn eff r m -^+m eff ^x(t) = F ext (t) , (3) 

which has the standard solution for displacement: 


We can relate this displacement to a physically observ¬ 
able quantity, i.e. the voltage produced at the mixer 
output of the microwave readout system: 

**(")= sxM iw) (£) ■ (5> 

We wish to solve this equation for the case of F ext origi¬ 
nating from applied RPF. The expression for this can be 
found by differentiating the standard interaction Hamil¬ 
tonian for an optomechanical system with respect to 
position 21 : 

Hint = - hg 0 a)a (b + 6 f ) , ( 6 ) 

where a) and a are the photon raising and lowering 
operators respectively, F and b are the phonon raising 
and lowering operators respectively, and go is the single 
photon optomechanical coupling: 

go = Gx zpf = 3 —^zpf, ^zpf = \ - - 7 ^— • (V 

dx V 2 

We can define the position operator x as 


Sx(cj) = 


1 F e xt fa ) 

-iujTm 


( 4 ) 


Here, m e ff is the effective mass of the bar in its resonant 
mechanical mode of frequency f2 m , T m = is the 

mechanical decay rate, and F ext is some applied force. 


x = xzpf (b + b^ , 
which will allow us to solve for F ext : 


( 8 ) 


Fext(t ) 


dHint 

dx 


KG a) a = h^-n cav 
dx 


(9) 
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Here, n cav is the number of photons inside the cavity, 
and can be solved by dividing the total energy inside the 
cavity 22 by the energy of a single photon, as follows: 


^cav 


-Ewgm 


E t 


singlephoton 


-E’si 


'singlephoton 


= huj 


-E'WGM — Fir 


Qe 


4/3i 


w (I + /81 +fo) 


1 + 


«32 (* 3 *)' 


( 10 ) 

where P^ nc is the power of the incident microwave sig¬ 
nal at frequency cj e , Q e the quality factor of the WGM, 
and pi and /?2 are the couplings between the WGM res¬ 
onance and the input and output microwave probes, re¬ 
spectively. 

Assuming we are driving the microwave input at the 
WGM resonance (cu = uj e ), we finally arrive at an expres¬ 
sion for the radiation pressure force inside the sapphire: 


Fext(t) | 


U)=U) e 


diU Pjnc(t)Qe 4/?i 
dx (1 + /?! + /3 2 ) 2 


( 11 ) 


As a side note, this equation is identical to that ob¬ 
tained by Locke et al . 23 when modelling the optomechan¬ 
ical system as an LCR circuit with a modulated capaci¬ 
tance. 

Assuming the incident power is modulated at fi m , we 
can derive an equation for the voltage output by the mi¬ 
crowave readout system: 


Su(Q m ) = ttV2tt 




2 

XWGM XmechSP 


( 12 ) 


where Xwgm = gf and X mech = 

From Eq. 12, we can observe that by measuring 
du/dP , one can determine a value for df /dx experimen¬ 
tally. The results of this measurement are shown in Fig. 
6 , and every other value in Eq. 12 can be directly mea¬ 
sured from the system. From the gradient of these two 
curves, df/dx 127 = 0.19 MHz//xm and df /dx 95 = 0.05 
MHz//im. 

The value for the 127 kHz mode is in very good agree¬ 
ment with the predicted value 6 % error), while the 
95 kHz mode is measured to be 40% smaller than its es¬ 
timated value (see Table I). This can be attributed to 
the assumption that there was no strain induced con¬ 
tribution to the 95 kHz transduction in the simulation. 
As the WGM used has ~ 95% of its E-field in the z- 
direction, any strain contribution would arise mainly in 
this axis, which would subtract from the value obtained 
from purely dimensional changes (see Eq. 1). For ex¬ 
ample, a strain contribution 1 / 10 th the magnitude of the 
127 kHz mode would result in the predicted value being 
in good agreement with the measured value. This arises 
in practice because the strain curve of the 95 kHz mode 



FIG. 7: Measured mechanical Q-factors as length of Nb 
suspension wire is altered. Lengths are measured from 
the angle made by the suspension. A speed of sound in 
Nb of 2200 m/s produces the best fits, in agreement 
with quoted values for the speed of transverse waves in 
a thin Nb rod at 20 °C. 


is not a pure odd function, and the WGM is not a pure 
even function about Z — 0 (see Fig. 3), as was assumed 
when treating the interaction of the two modes with Eq. 
1 , which, in addition, is only an approximation to first 
order. Some small perturbations of these two symme¬ 
tries will result in a small strain contribution to df/dx 95, 
which will subtract from the purely dimensional contri¬ 
bution. 


IV. OPTIMISING MECHANICAL Q-FACTORS 

Energy dissipated through the suspension is generally the 
dominant loss mechanism in optomechanical systems, as 
long as they are held under vacuum. 

The wire-loop suspension used to hold the SB is a 
method traditionally used by the gravitational wave com¬ 
munity, which has been shown to achieve the highest me¬ 
chanical Q-factors over other suspension schemes 24 (see 
the comparison of suspension techniques depicted in Fig. 
9). Braginsky et a /. 24 modelled the losses in such a sus¬ 
pension as; 


1 

Q~ s 


‘ZpwAwlw X ext x 

'WleffQw x q 



^mfiL 

x wQvl 


2 

— COS 


2L / ,cj r7 



(13) 


where p w , A WJ l w , Q w and v w are the density, cross- 
sectional area, length, mechanical Q-factor and velocity 
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FIG. 8: Mechanical Q-factor of the 127 kHz and 95 kHz 
modes vs pressure. Residual gas damping limits as 
predicted by equation (15). 


of sound of the wire, respectively, and xo and x ext 
are the amplitude of vibration of the resonator and 
the amplitude at the point of contact with the wire, 
respectively. Equation 13 predicts a periodic variation in 
Q-factors as the length of the wire (shown in blue in Fig. 
4) is changed. However, Eq. 13 has only been shown to 
provide qualitative agreement with experimental results, 
as the loss mechanism is far more complicated than this 
simple model. 

Two different wires were trialled for the suspension sys¬ 
tem; an 80 /i m diameter tungsten wire and a 125 p m 
niobium wire, with better results being achieved with the 
latter due to its higher intrinsic Q w . The niobium clamp 
(shown in Fig. 4) is necessary to keep enough friction 
around the neck of the SB to maintain balance, but can 
also be used to alter the length of the wire by moving its 
position. The clamps are made by cutting small lengths 
of 0.65 mm diameter niobium tubing. Tests were made 
with both clamped and unclamped pieces of tubing, with 
the former producing higher Q-factors due to the bound¬ 
ary condition in the wire resonance being fixed. 

Results of varying wire length are presented in Fig. 7. 
Each experimental data point represents an entirely new 
wire suspension, as the clamping technique used fixes the 
wire’s length. As such, it is possible for surface conditions 
and the loop’s position along the z-axis to change be¬ 
tween measurements. Therefore each of the data points 
in Fig. 7 represents a lower limit of Q-factor for any given 
wire length. 

Braginsky et al. 24 states that there is only a qualita¬ 
tive agreement between Eq. 13 and the experimental 


data, with measured Q-factors being larger than the pre¬ 
dicted values by almost a full order of magnitude. This is 
mainly due to resonator vibrations at the point of contact 
with the suspension wire being only partially transmit¬ 
ted to the wire. As such, we insert an additional term 
inside the square brackets of Eq. 13 to represent cou¬ 
pling between the wire and crystal. It takes the form of 
l/£. It should also be noted that the second term inside 
the square brackets of Eq. 13 is much less than 1 and 
can therefore be neglected. The experimental results are 
then fitted with the following equation: 


Qs = 


TH'eiiQ'U 


^Pw^-w^w %ext 


1 (2l w UJ rn 


1 + - — COS 


, ( 14 ) 


which is shown in black in Fig. 7 and provides a good 
bound of the measured Q-factors given the fitted values 
of Q w = 8.5 x 10 5 and f = 0.31. 


Fig. 8 shows the dependence of mechanical Q-factors on 
pressure, as predicted by the following equation: 25 

1 ~ PA 
Qgas ATciJo 

where P is pressure (in Pascals), A the surface area 
of the mechanical resonator, M its mass and p is the 
mass of molecule comprising majority of the gas in 
question. As can be seen from Fig. 8, below a certain 
threshold pressure, residual gas damping is no longer 
the dominant loss mechanism, with Q-factors of close to 
10 8 demonstrated. 



Throughout the literature regarding wire loop suspen¬ 
sions, it is often noted that by applying a small layer of 
lubrication (most commonly animal fat is used) to the 
surface of the mechanical resonator at the point of sus¬ 
pension, significant improvements in Q-factors can be 
achieved 13,24,26 . Previously, at room temperature, Q- 
factors on the order of 10 8 and Q m x / > 6 x 10 12 have 
only been achieved with some form of lubrication, and 
this is therefore a strategy that could be implemented 
in the immediate future of this work. However even 
without lubrication, the present work corresponds to a 
Qm x / product of 10 13 , equivalent to the best mea¬ 
sured values of sapphire acoustic systems, as depicted 
in Fig. 9. This product is considered a figure of merit 
for acoustic resonant systems, as it is a direct measure 
of the degree of decoupling from the thermal environ¬ 
ment. Specifically, Q m x / > k^T/h is the condition 
for neglecting thermal decoherence over one mechanical 
period, and is considered the minimum requirement for 
quantum optomechanics 21 . At room temperature, this 
lower limit is equal to 6 x 10 12 , and is represented in Fig. 
9 by the red dashed line. 

Importantly, our results are consistent with the 
best reported values measured at room temperature 
in sapphire, most likely limited by the quality of the 
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crystal 13,24 . Thus, this result demonstrates that equiv¬ 
alent performance to state-of-the-art sapphire acoustic 
systems can be achieved with a modified cylindrical 
structure. All other results presented in Fig. 9 are 
achieved using modes in right cylinders as shown in the 
right-hand side limit of Fig. 2. This is the first time 
optomechanics in a sapphire bar has been attempted 
with a modified architecture. 
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FIG. 9: Reported mechanical quality factors and 
resonant frequencies of published experiments with 
sapphire bar optomechanical systems at room 
temperature: 1. (Rowan et al , 2000) 26 , 2. (Locke et a/., 
2001) 13 , 3. (Braginsky et a/., 1985) 24 , 4. (Numata et 
al , 2000) 27 , 5. (Uchiyama et al , 1999) 28 . The 
experiments are differentiated by their suspension 
techniques. The black dashed line represents the Q x / 
product of the SB. 

By cooling the SB, the coupling between the suspen¬ 
sion and crystal will decrease, resulting in the suspension 
losses decreasing. Q-factors on the order of 10 9 have been 
achieved using wire suspended sapphire at liquid helium 
temperatures 24 . So whilst a maximum value of ^ 8 x 10 7 
has been achieved thus far, we are optimistic of improv¬ 
ing this. 

V. MICROWAVE READOUT SYSTEM 

A. Detection of Frequency Fluctuations 

As described in the previous sections, a freely sus¬ 
pended sapphire bar is a unique physical object capa¬ 
ble of supporting very high quality resonances, both me¬ 
chanical and electromagnetic. It offers an opportunity 
to conduct a detailed study of the bar’s mechanical res¬ 
onances via their influence on the microwave resonances. 


In addition, one can investigate the inverse effects of the 
microwave readout on the SB’s mechanical properties in¬ 
cluding the degeneration (cold damping) and regenera¬ 
tion (parametric excitation) of elastic vibrations 24,29 . 

This section describes a microwave readout system for 
monitoring the vibration state of a freely suspended sap¬ 
phire bar. We expected that the sensitivity of the mi¬ 
crowave readout would be sufficiently high to permit the 
very first observations of the bar’s mechanical resonances 
excited by thermal fluctuations of the sapphire crystalline 
lattice. 

There are at least two techniques, which can be em¬ 
ployed for the microwave assisted detection of elas¬ 
tic vibrations of the SB. Both techniques are widely 
used in the field of oscillator frequency stabilisa¬ 
tion for high-resolution measurements of fast frequency 
fluctuations 30,31 . In one case, the SB can be configured 
as a dispersive element of a microwave frequency discrim¬ 
inator. When driven from a fixed frequency signal source, 
the SB acts as a frequency-to-voltage converter producing 
voltage varying synchronously with resonant frequency of 
given WG-mode excited by the source. As a result, the 
task of analysing spectrum of a microwave signal is re¬ 
duced to computing the Fast Fourier Transform of a sam¬ 
pled voltage at the discriminator output. Alternatively, 
the SB can serve as a band-pass filter of a self-exciting 
microwave loop oscillator 32 . This would “imprint” the 
spectrum of the SB mechanical normal modes on to the 
spectrum of the microwave signal. Once again, a fre¬ 
quency discriminator can be used to convert frequency 
fluctuations of the microwave signal into synchronous 
fluctuations of voltage. Performance-wise, both of the 
above-mentioned techniques are identical; in each case 
the useful frequency fluctuations associated with elastic 
vibrations of the SB must compete with the same spu¬ 
rious fluctuations of the microwave readout electronics; 
and each technique, given some optimisation, could po¬ 
tentially lead to spectral resolution close to the Standard 
Thermal Noise Limit 33 . 


B. Microwave Sapphire Bar Oscillator 

Figure 10 shows a schematic diagram of a microwave 
loop oscillator based on the SB. The oscillator is 
frequency-locked to a given WGM of the SB to facilitate 
detection of the useful frequency fluctuations. The spec¬ 
trum of the free-running oscillator at Fourier frequencies 
of interest (^ 100 kHz) is completely dominated by 1 // 
phase noise of the microwave loop amplifier. These spu¬ 
rious fluctuations are suppressed by the frequency lock 
(or frequency control) system as explained below. 

The key element of the oscillator frequency control sys¬ 
tem is an ultra-sensitive frequency discriminator. It con¬ 
sists of a microwave Mach-Zehnder interferometer with 
the suspended sapphire bar resonator and a phase sensi¬ 
tive readout system featuring low-noise amplifier (LNA) 
and a double-balanced mixer (DBM1 in Fig. 10). The 








Filtered & 
low-power <—« 
microwave 
output 

Sine wave at 90... 150 kHz 
to excite mechanical 
resonances in 

sapphire bar VVV >{ CL 
VCA 


Amplifier 


To external 
readout 
system 


Legend: 

BPF - band-pass filter 

CRL - circulator 

DBM-double-balanced mixer 

DC - directional coupler 

F - loop filter 

HC - 3 dB hybrid coupler 

LNA - low-noise amplifier 

VCA - voltage-controlled attenuator 

VCP - voltage-controlled phaseshifter 

cp - mechanical phaseshifter 


DC 



Modulation 


FIG. 10: Schematic diagram of a frequency stabilized microwave oscillator based on suspended sapphire SB 
resonator excited in electromagnetic modes of whispering gallery. 


microwave signal reflected from the SB interferes destruc¬ 
tively with a fraction of the incident signal at the inter¬ 
ferometer “dark port”. This cancels the carrier of the 
difference signal while preserving noise modulation side¬ 
bands resulting from 1 //-noise of the microwave loop am¬ 
plifier. The residual noise at the dark port is amplified 
and demodulated to the baseband, producing an error 
voltage proportional to oscillator frequency fluctuations. 
The error voltage, after appropriate filtering, is applied 
to the electronic phase-shifter (VCP) in the microwave 
loop. This steers the oscillator frequency to that of the 
resonator, or more precisely, to the frequency at which 
the carrier was suppressed, which is typically well within 
the resonator bandwidth. The frequency discriminator, 
loop filter and VCP form the frequency control loop. If 
the control loop gain is sufficiently high, the fidelity with 
which oscillator frequency follows that of the resonator 
is determined only by technical fluctuations in the elec¬ 
tronics of the frequency discriminator. In this respect, 
interferometric frequency discriminators are far superior 
to their conventional counterparts, as they exhibit effec¬ 
tive noise temperature close to the ambient temperature, 
and are capable of handling much higher power levels 33 . 

It is not difficult to understand why frequency discrim¬ 
inators with the highest sensitivity are based on spindle- 
shaped sapphire crystals fixed rigidly inside the protec¬ 
tive metal shields. The main reason for this is the low 
vibration sensitivity of such resonators, which makes the 


task of carrier suppression fairly straightforward. The 
situation is different in the case of the SB whose rock¬ 
ing motion upsets both amplitude and phase balance of 
the interferometer. To cope with the SB rocking motion, 
the oscillator in Fig. 10 features an additional feedback 
control system charged with the task of minimising am¬ 
plitude mismatch between two signals interfering at the 
dark port. The error voltage for this feedback system is 
produced by the second mixer (DBM2) tuned in quadra¬ 
ture relative to the mixer of the frequency control loop. 
The voltage-controlled attenuator (VCA) in the interfer¬ 
ometer arm completes the feedback loop acting as an 
actuator of the amplitude mismatch control system. 

The quadrature tuning of both DBMs is required to 
avoid the cross-talk between the frequency and amplitude 
control loops. Yet, the tolerances of such tuning proved 
to be not very stringent. For each control system, it 
was sufficient to roughly maximise the amplitude of the 
error signal (in response to some deliberately introduced 
perturbation) before closing the feedback loop. 

An additional VCA in Fig. 10 is used to modulate 
the amplitude of microwave signal. When the modula¬ 
tion frequency coincides with the frequency of mechanical 
resonance, parametric excitation of the normal mode en¬ 
sues. This VCA was introduced when the refined suspen¬ 
sion system made piezoelectric excitation of mechanical 
resonances ineffective. 
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Legend: 

BPF - band-pass filter 

CRL - circulator 
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FIG. 11: Schematic diagram of interferometric frequency discriminator with automatic carrier suppression. 


C. Inteferometric Frequency Discriminator 

Approximately half of the power generated by the SB 
oscillator (Fig. 10) is diverted to the external frequency 
discriminator shown schematically in Fig. 11. It is prac¬ 
tically identical to the previously described built-in dis¬ 
criminator of the SB oscillator - it also contains a mi¬ 
crowave interferometer with actively controlled balance 
to cope with frequency variations of the input signal in¬ 
duced by the rocking motion of the SB and changes of 
ambient temperature. 

The dispersive element of the external frequency dis¬ 
criminator is a cylindrical metal cavity with a Q-factor 
of 15,000 and frequency tuning range of approximately 2 
GHz. The wide tuning range of the cavity enables easy 
switching from one WGM of the sapphire bar to another, 
if one needs to investigate how a WGM of different po¬ 
larisation and azimuthal number responds to vibration. 

Fig. 12 shows spectra of phase fluctuations of various 
signal sources measured with the external frequency dis¬ 
criminator. The measurements were made at 9.774 GHz 
corresponding to the excitation of the WGHi 5 ? i ? i mode 
of the SB. 

The power spectral density of phase fluctuations 
was inferred from that of the voltage noise S u via the 
following relationship: 


t? n 




(16) 


where F is the Fourier frequency, A/z, is the half- 
loaded bandwidth of the cavity resonator and Kfd is 
the frequency discriminator “DC sensitivity” measured 
at F « Afh- 

The bottom trace in Fig. 12 shows the fit to the noise 
floor of the frequency discriminator expressed in the sin¬ 
gle sideband (SSB) units of dBc/Hz. The noise floor 


was measured with the cavity resonator replaced by a 
50 D termination. Next, we measured phase noise of a 
commercial frequency synthesizer (Agilent E8257C). The 
idea was to verify that the voltage-to-phase conversion 
procedure we followed was correct (our results proved to 
be consistent with the specs of the Agilent instrument). 
Finally, we characterised the phase noise of the SB os¬ 
cillator, both in the free-running and frequency-locked 
regimes. At Fourier frequencies 50 kHz < F < 500 kHz, 
we observed more than 20 dB of phase noise suppression 
when the frequency lock was engaged. 

As follows from Fig. 12, the SSB power spectral den¬ 
sity of spurious phase fluctuations is -165 dBc/Hz at 
F=100 kHz. The corresponding level of rms frequency 
fluctuations is 


5f = Fy/S^F) - 8 x 10 - A Hz/yW~z (17) 

Recalling the frequency-displacement sensitivity of the 
sapphire bar ( f z ~ 0.18 MHz/gm) yields the displace¬ 
ment noise floor: Sz ~ fif /f z ~ 4 x 10“ 15 m/V^z. This 
corresponds to the highest displacement sensitivity ever 
reported in experiments with sapphire bar resonators as 
displacement transducers. Yet, the level achieved is still 
approximately twice as large as the 127 kHz mode’s elas¬ 
tic vibrations driven by thermal noise, hence Brownian 
motion is currently unobservable with the present con¬ 
figuration. 

At this stage, two questions can be posed: “What lim¬ 
its spectral resolution of the current experimental setup 
and how can it be improved?” In answering these ques¬ 
tions, we, first of all, can single out two factors, which 
are almost equally responsible for the present level of dis¬ 
placement sensitivity. One of them is the relatively high 
level of phase noise of the frequency-locked microwave 
oscillator at F ^ 100 kHz. This is because the current 
frequency control system starts loosing gain at Fourier 
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FIG. 12: SSB phase noise spectra of various microwave 
signal sources measured with interferometric frequency 
discriminator with automatic carrier suppression. 
Power of the input signal is 19 dBm. 


frequencies above 10 kHz. Solving this problem will in¬ 
volve the design of a new low-pass filter of the frequency 
feedback loop based on operational amplifiers with gain- 
bandwidth product exceeding 1 GHz. 

The second limitation arises from the comparatively 
low electrical Q-factor of the hollow metal cavity. We 
believe, that at least an order of magnitude improvement 
in displacement sensitivity can be gained by replacing 
the hollow metal cavity with a sapphire loaded cavity 
resonator. 

One “side effect” of resonator substitution is the loss 
of the wide frequency tunability. Yet, some residual 
tunability would remain owing to the relatively high 
sensitivity of sapphire resonators to temperature (df/dT 
~ 0.5-0.7 MHz/K depending on the type of WGM 
used). Another complication of resonator substitution is 
related to the narrow bandwidths of sapphire resonators. 
An additional control system would be required to keep 
the sapphire resonator “in sync” with the incoming 
signal. We plan to address this issue by controlling 
microwave power dissipated in the sapphire crystal 
as in 34,35 , where such a technique was used to enable 
the phase referencing of a “slave” sapphire oscillator 
to the “master”. It should be remembered that the 
measurement sensitivity improves as a square root of 
the power of the SB oscillator. 

A major advantage of the SB’s dumbbell-shaped archi¬ 
tecture is that it is in fact two electromagnetic oscillators, 
which are both undergoing the same mechanical fluctua¬ 
tions. One can easily imagine a cross correlation scheme 


in which two readout systems as described above were 
constructed around either end of the SB. The primary 
mixer outputs of these readouts could then be cross corre¬ 
lated to eliminate uncorrelated electronic noise and boost 
the correlated mechanical signal, which would potentially 
allow never-before-seen levels of displacement sensitivity 
in such a system. 


VI. PARAMETRIC EFFECTS 

Parametric behaviour will always arise when a driven 
resonant system is coupled to a second resonant system. 
A phase difference between the mechanical oscillator 
and the microwave resonance will produce either a 
damped or driven system. As such, one can expect 
changes in the mechanical quality factor of a resonant 
optomechanical system as the “pump” signal is detuned 
from resonance. These effect have been well described in 
the past with various types of transducers 17,21,23,36,37 . 
The phenomenon, referred to as the “optical spring” 
effect, results from a Stokes/anti-Stokes process, 
whereby pumping above the optical resonance produces 
additional phonons; increasing the mechanical mode’s ef¬ 
fective temperature and Q-factor, whilst pumping below 
resonance removes phonons from the system, effectively 
increasing losses and cooling the system 38 . There should 
also be a corresponding mechanical frequency shift 
associated with microwave pump detuning, however it 
is predicted that it would be on the order of mHz, far 
below the resolution of the measurement technique used 
here, which relies upon a vector signal analyser’s “peak 
trace” function to record amplitude vs. time. 

Figure 13(a) demonstrates the “optical spring” 
effect for the 127 kHz mode of the SB resonator. These 
measurements were taken using a method depicted in 
Fig. 13b, in which whilst the mechanical resonance rings 
down, the microwave detuning is changed from zero, 
and the resulting change in time constant is measured. 
Positive detuning results were also taken and the trials 
repeated. The transient response that can be seen 
immediately after a frequency shift in Fig. 13b is a 
result of the readout system requiring a re-balance. 

Within the explored detuning range, A co e < £2 m , 
hence one expects an approximately linear relationship 
between Q-factor and detuning. The reason for this 
limited range is due to the restrictions placed by the 
bandwidth of the frequency control system acting on the 
loop oscillator (see Section V). The loop oscillator and 
control system used to readout the mechanics of the SB 
complicate the microwave system beyond a simple LCR 
model. 

The presence of parametric effects is a promising result 
for the system, as it is only through resolved sideband 
cooling 38 that a system such as the SB resonator could 
overcome thermal noise to reach a quantum limited 
state 21,38 . 
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To date, we have achieved excellent agreement be¬ 
tween modelled and measured frequencies of a variety of 
mechanical modes of the SB resonator, and developed 
a robust method of mechanically exciting it via a 
piezoelectric shaker and via radiation pressure. The 
later method has also provided a method of calibrating 
the transducer’s frequency sensitivity to displacement, 
df/dx , which has been measured to good agreement 
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FIG. 13: (a) Q-factor vs input microwave frequency 
detuning. The circled data points are obtained from 
analysing (b); a real time ringdown measurement of the 
127 kHz mechanical mode as negative detuning is 
increased during the measurement. The natural log of 
the amplitude is taken such that the time constant r 
will be equal to the inverse of the gradient of a linear fit 
for each detuning value’s data. 


with predicted results. The mechanical Q-factors have 
been optimised to a point at which they agree with 
maximum reported values at room temperature within 
limitations set by the quality of the sapphire, as does 
the Q m x / product. Parametric back-action resulting 
in mechanical damping and excitation have also been 
observed. Finally, a novel microwave readout system has 
been constructed, which provides extremely low phase 
noise performance. This work provides an important 
enabling step for the next generation of kilogram scale 
mechanical oscillator experiments designed to measure 
and test the standard quantum limit and potentially 
investigate the nature of quantum gravity, thus proving 
that novel architectures of sapphire systems can still 
produce state-of-the-art results. 
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